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We discuss production of heavy partciles during reheating. We find that the very 
energetic inflaton decay products can contribute to the production of massive stable 
particles, either through collisions with the thermal plasma, or through collisions with 
each other. If such reactions exist, the same massive particles can also be produced 
directly in inflaton decay, once higher-order processes are included. We show that these 
new, non-thermal production mechanisms often significantly strengthen constraints on 
the parameters of models containing massive stable particles. 



After inflation W, coherent oscillations of the inflaton dominate the energy density of the Universe. 
At some later time these coherent oscillations decay to the fields to which they are coupled, and their 
energy density is transferred to relativistic particles; this reheating stage results in a radiation-dominated 
Friedmann-Robertson-Walker (FRW) universe. 

Until a few years ago, reheating was treated as the perturbative, one particle decay of the inflaton with 
decay rate T^, leading to the simple estimate Tr ~ iXdMpundi) 1 ^ 2 for the reheat temperature where 
Mpianck = 2.4 x 10 18 GeV represents the reduced Planck mass. It has been noticed in recent years that 
the initial stage of inflaton decay might occur through a complicated and non-perturbative process called 
parametric resonance ||. However, it is generally believed that an epoch of (perturbative) reheating 
from the decay of massive particles (or coherent field oscillations, which amounts to the same thing) is 
an essential ingredient of any potentially realistic cosmological model Q . In what follows we generically 
call the decaying particle the "inflaton" . However, it should be clear that our results are equally well 
applicable to any other particle whose (late) decay results in entropy production. 

Even before all inflatons decay, the decay products form a plasma which has the instantaneous tem- 
perature 



where H is the Hubble parameter and denotes the number of relativistic degrees of freedom in the 
plasma |J. This temperature reaches its maximum T max soon after the inflaton field starts to oscillate. 
In addition to this thcrmalizcd plasma there are inflaton decay products with energy ~ which will 

eventually come into equilibrium with the thermal bath. 

In recent years several mechanisms have been put forward for creating very heavy, even superheavy, 
particles in cosmologically interesting abundances [^j. Here we will focus on production of very massive 
particles from various processes, including a thermal bath, during perturbative reheating. Note that 
particle production from other sources, if present, would further strengthen the bounds which we will 
derive as they simply add to production from mechanisms discussed here. 



In Rcf. out of equilibrium production of x from scatterings of "soft" particles in the thermal bath 
(with energy E ~ T) is studied and the final result is found to be (the superscript "ss" stands for \ 
production from "soft-soft" scattering) 
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(x not in equilibrium). (2) 

In the opposite situation, with \ being initially in equilibrium with the thermal bath, today's x relic 
density is given || by 
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where the exponent a — (1) if XX annihilation proceeds from an S (P) wave initial state. The freeze-out 
temperature is now given by Xf = m x /T{ ~ log^.OSg.T^a! 2 £f '^"-MpianckTR/m 3 ). 

However, as mentioned earlier, "hard" particles with energy E ~ 3> T are continuously created 

by inflaton decay for H > These particles eventually thermalize with the bath, but this takes a 
finite amount of time. The presence of hard inflaton decay products can therefore affect heavy particle 
production in two ways. Firstly, xs can be produced from 2^2 scatterings of a hard particle off cither 
soft particles in the thermal bath (if kinematically allowed) , or off other hard particles Moreover, xs 
might be directly produced from inflaton decay |10| . 

A. Particle production from hard— soft scatterings 

In order to estimate the rate of heavy particle production from the "hard" inflaton decay products, 
we also have to know the time needed to reduce their energy from a value ~ to a value near T. 

As shown in Ref. [jnj, 2^2 scattering reactions are not very efficient in this respect. The reaction 
rate is large, but the average energy loss per scattering is only 0(T 2 /m^), giving a slow-down time of 

order [a 2 T 2 /mA 1 (up to logarithmic factors). On the other hand, inelastic 2 — > 3 reactions allow large 
energy losses (in nearly collinear particles) even if all virtual particles only have virtuality of order T. 
The slow-down rate is thus given by: 

r,,^ **r (&)"*. (4) 

Next let us estimate the rate for x pair production from hard-soft scatterings. This process is kine- 
matically allowed so long as ET > 4m 2 ,, where E is the energy of the hard particle so that the square of 
the center-of-mass energy is typically a few times ET. The hard particle initially has energy E ~ m^/2 
and average number density ~ g*T A / (3m^,), just after its production from inflaton decay. On the 
other hand, the rate for x production from hard-soft scatterings is approximately given by 

, / a 2 aa? \ 
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The two contributions in (|^) describe 2 — > 2 reactions with squared center-of-mass energy ~ m^T and 
"radiative return" 2 — > 3 reactions, respectively; in the latter case the hard particle emits a collinear 
particle prior to the collision, thereby reducing the effective cms energy of the collision to a value near 
m x . In order to make a safe (under)estimate we choose the temperature To = 2T t hr for presenting our 
results; note that the \ pair production cross section at threshold, s = 4m 2 , is suppressed kinematically. 
Our final results for the contribution of hard-soft collisions to the x re hc density are M: for T < Tr: 
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in the opposite situation, we have 
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(T > T R ) (7) 



B. Particle production from hard— hard scatterings 



If To > 7r, we should also consider x production from scattering of two hard particles. Collisions 
of these particles with each other can produce x pairs if m x < m$/2. Note that this constraint is 
independent of the temperature. On the other hand, it's also possible that To > T max , in which case 
hard-soft scattering (and soft-soft scattering) does not produce any \ particles. 

The rate of x production from hard-hard scattering is quadratic in the density of hard particles. 
Therefore we can not use our earlier approximation of the density of hard particles produced in one 
Hubble time in the presence of a thermalized plasma, since the actual density n^t) at any given time 
will be much smaller than this. In a plasma with temperature T a hard particle will only survive for a 
time ~ l/T slow , see eq.([|). For T max > T > T R , the production of hard particles from inflaton decays 
and their slow-down will be in equilibrium, i.e. the instantaneous density nh(t) — 2Tdn < j > (t)/T a i ovr , where 
is the density of inflatons. By taking into account x production prior to the bulid up of a thermalized 
plasma, our final estimate will be 

nfh 2 ~ 6 ■ io 27 ■ ( 1/2 o xx . (8) 

In Fig. (2) we present three numerical examples to compare the significance of hard-soft and hard-hard 
scatterings with that of soft-soft scatterings. A detailed discussion on this figure and features observed 
in it can be found in Ref. 0. 



C. Particle production from inflaton decay 

We now discuss the direct production of x particles in inflaton decay whose importance has recently 
been noticed ]lO| j. Let us denote the average number of x particles produced in each <fi decay by B(cj) — > x). 
The x density from <fi decay can then be estimated as: 

ftdecay^ „ 2 . 1O 8 B(0 _> (9) 

* rritj) 1 GeV 

Eq. (||) holds if the % annihilation rate is smaller than the Hubble expansion rate at T ~ Tr. 

We now discuss estimates of B((f> — > x)- This quantity is obviously model dependent, so we have 
to investigate several scenarios. The first, important special case is where x is the LSP. If is large 
compared to typical visible-sector superparticle masses, <f> will decay into particles and superparticles 
with approximately equal probability |10[. Moreover, all superparticles will quickly decay into the LSP 
and some standard particle(s). As a result, if x is the LSP, then B(<fi — > \) — lj independently of the 
nature of the LSP. 

Another possibility is that the inflaton couples to all particles with more or less equal strength, e.g. 
through non-renormalizable interactions. In that case one expects B(cf> — > x) ~ ~ 1/200. However, 
even if <f> has no direct couplings to x, the rate (Q) can be large. The key observation is that x can be 
produced in <f> decays that occur in higher order in perturbation theory whenever x can be produced from 
annihilation of particles in the thermal plasma. In most realistic cases, cf> — > ffxx decays will be possible 
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FIG. 1: Sample diagram for x production in four-body inflaton decay. 

if x has gauge interactions, where / stands for some gauge non-singlet with tree-level coupling to 0. A 
diagram contributing to this decay is shown in Fig. (1). Note that the part of the diagram describing 
XX production is identical to the diagram describing XX // transitions. This leads to the following 
estimate: 

Cioi ( 4ml \ 2 ( 2m v \ ' 
967r d y m i J \ m 4> ) 

where C4 is a multiplicity (color) factor. The phase space factors have been written in a fashion that 
reproduces the correct behavior for m x — » m^/2 as well as for m x — > 0. This estimate provides a lower 
bound on B(cj) — > x) under the conditions assumed for our calculation of fi!j s and fi^h. whenever a primary 
inflaton decay product can interact with a particle in the thermal plasma, or with another primary decay 
product, to produce a XX pair, cf> — > x four-body decays must exist. 

Occasionally one has to go to even higher order in perturbation theory to produce x particles from (f> 
decays. For example, if x has only strong interactions but <j> only couples to SU (3) singlets, XX pairs can 
only be produced in six body final states, — > ffqqxX- A representative diagram can be obtained from 
the one shown in Fig. (1) by replacing the x lines by quark lines, attaching an additional virtual gluon 
to one of the quarks which finally splits into XX- The branching ratio for such six body decays can be 
estimated as 

C 6 ala 2 I . 4ml \ I '. 2m, 



Finally, in supergravity models with explicit (supersymmetric) x niass term there in general exists a 
coupling between <f> and either x itself or, for fermionic x> to its scalar superpartner, resulting in the 
estimate Q 

167r \/^ M Planck T R V m 4> J 

where v denotes the vacuum expectation value of the inflaton at the true minimum of its potential. 



III. DISCUSSION 

The production of x particles from inflaton decay will be important for large m x and large ratio m x /T^, 
but tends to become less relevant for large ratio m < j > /m x . Even if m x < T max , x production from the 
thermal plasma (|^) will be subdominant if 

B{<t> ^ X) > f IOOTr \ 6 m$ 1 TeV 
a\ \ m x ) m x m x ' 
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Note that the first factor on the r.h.s. of ( |l3| ) must be < 1CP 6 in order to avoid over-production of \ 
from thermal sources alone. 

In |hJ we showed that the decay contribution (^|) by itself leads to very stringent constraints on 
models with massive stable \ particles. In particular, charged stable particles with mass below ~ 100 
TeV seem to be excluded, unless m x > In case of a (neutral) LSP with mass around 200 GeV, 

the overclosure constraint implies tti^/Tr > 4 • 10 10 , i.e. a very low reheat temperature, unless \ wa s 
in thermal equilibrium below Tr; recall that B(<f> — > \) — 1 m this case. Finally, if ~ 10 13 GeV a 
"wimpzilla" with mass m x ~ 10 12 GeV will be a good Dark Matter candidate only if it has a very low 
branching ratio, B((f> — > \) ~ 5 • 10~ 8 GeV/Za, i.e. if its couplings to ordinary matter are very small. 

Many of the results presented here are only semi-quantitative. Unfortunately in most cases significant 
improvements can only be made at great effort. For example, a proper treatment of the slow-down of 
primary inflaton decay products would require a careful treatment of the full momentum dependence of 
the particle distribution functions. On the other hand, our estimates of B(cf> — > x) should be quite reliable 
if m x > Tr (which is required for x n ot to have been in thermal equilibrium at Tr); even for many-body 
decays, details of the matrix elements should change our estimates only be 0(1) factors. Fortunately this 
is often also the most important of the new mechanisms for the production of massive particles at the 
end of inflation. 
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FIG. 2. The relic density \ particles would currently have if they are absolutely stable is shown as a function 
of m x for couplings a = 0.05, a x = 0.01, and a) (T R ,m ) = (10 8 GeV, 10 13 GeV), b)(10 5 GeV, 10 13 GeV), 
c)(3 MeV, 10 s GeV). The soft-soft, hard-soft and hard-hard contributions are shown by the dotted, short dashed 
and long dashed curves, respectively, while the solid curves show the sum of all three contributions. 
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